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Геоинформационная система анализа 
климатических и физико-географических 
параметров для исследования ареалов 
лесных пород 


В статье рассматривается геоинформационная система анализа климатических (средняя температура) 
и физико-географических (абсолютная высота) параметров для исследования ареалов лесных пород. 
Предложен подход к анализу параметров, который позволяет моделировать сценарии изменения 
климата и прогнозировать влияние изменения климата на ареалы лесов. 


Введение 


Изменение климата является наиболее важной и сложной проблемой в сфере 
охраны окружающей среды, которая возникла перед человечеством за последний 
век. Подписание Рамковой Конвенции ООН об Изменении Климата представителями 
150 стран свидетельствует о том, что изменение климата существенно угрожает 
экологии Земли и экономическому развитию человечества [1]. Известно, что за 
последние сто лет наша атмосфера заметно потеплела и этот процесс продолжается. 
Климат планеты постоянно меняется. Изменение климата ведет к изменению ареалов 
всех видов растений, в том числе и к изменению территории распространения лесов [2]. 
Поэтому актуальной научной задачей является разработка математического, 
алгоритмического и программного инструментариев для анализа влияния климатических 
параметров на ареалы лесных пород. 

Целью данной работы является разработка математического аппарата и 
геоинформационной технологии для формирования базы климатических данных об 
ареалах лесных пород для западной и центральной части Украины. 


Формирование цифровых карт лесов, климатических 
и физико-географических параметров 


Для проведения анализа климатических и физико-географических параметров 
(температура, влажность, количество осадков, абсолютная высота и др.), какие присущи 
для ареалов лесов, использованы средства геоинформационной системы МарШ®Ю, которая 
дает возможность редактирования, отображения и анализа географической информации. 
При помощи такой геоинформационной системы можно не только создавать модельное 
изображение территории и разнообразных географических объектов, но и задавать 
произвольную информацию табличного типа для каждого из них [3]. Инструментарий 
геоинформационной системы является эффективным и удобным в проведении анализа 
климатических и физико-географических параметров для ареалов лесов. 
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Для проведения исследования необходимо сформировать карты ареалов лесных 
пород, климатических и физико-географических параметров и соответствующие 
таблицы с необходимыми данными. 

Анализ проведен для западных и центральных областей Украины (Винницкой, 
Волынской, Житомирской, Закарпатской, Ивано-Франковской, Киевской, Львовской, 
Ровенской, Тернопольской, Хмельницкой, Черкасской, Черновицкой), поскольку именно 
в этих областях наибольшая площадь ареалов лесов [4], [5]. Для осуществления расчетов 
выбраны такие физико-географические и климатические параметры, как абсолютная 
высота и средняя температура. 

Как входные данные используется цифровая карта потенциальной растительности, 
созданная с помощью средств МаршЮ. Под потенциальной имеется в виду 
растительность, которая была бы на соответствующей территории в настоящее 
время, если бы не было хозяйственной деятельности [6], [7]. Цифровая карта (рис. 1) 
отображает расположение растительности на территории анализируемого региона. 


[Е] Карпатские горные и предгорные леса (Ель, пихта, бук, дуб) 
= Широколистные леса (Бук, дуб, граб, клен, липа) 
Е Луговые степи (лаковые, разнотравье) 


Е | Хвойные и широколистные леса (Сосна, дуб, береза, граб, ольха, липа) 


Рисунок 1 — Карта потенциальной растительности для западных и центральных 
областей Украины 


Для формирования карты абсолютных высот региона использована цифровая карта 
высот Национального атласа Украины. Построение карты высот для анализируемого 
региона происходит по такой схеме: 

1) наложить на цифровую карту высот карту рассмотренных областей Украины, 
которая содержит информацию о границах областей и их другие параметры; 

2) карту высот расположить как верхний слой и установить возможность ее 
редактирования; 
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3) выделить все объекты карты высот, использовав опцию для выбора объекта, 
который изменяется; 

4) установить карту высот как нижний слой, выделить области Украины и удалить 
внешнюю часть карты высот, то есть ту часть, которая не входит в предварительно 
выбранные области. 

Такую процедуру можно сделать с помощью запроса языком МарВаз1с: 


Мар Егот ВеПеЁ 21, ВезВоипа_ НаОКг 

Зе Мар Гауег 1 ЕацаЫе Оп 

зеЁ тар гедга\ оп 

ЗеЕ Тагое! Оп 

зеё тар гедга\у о 

Зе Мар Огаег 2,1 

зеё тар гедга\ оп 

ОБесё$ ПиегзесЕ шо Тагоеё аа ПинаП=ПинаПо, ОТО=ОТО, ЗВееЙ ауег=ЗБееГ ауег, 
СодеТоро=СодеТоро, Нез|НКе!=НеюНВВе!|, Нез» АБ з=НеюВАЬ$ 


Полученная цифровая карта отображает расположение изолиний абсолютных высот 
на рассмотренной территории Украины (фрагмент карты с изолиниями показан на рис. 2). 


— 150 — Изолинии абсолютных высот 


Рисунок 2 — Фрагмент карты изолиний абсолютных высот 


Поскольку на большей части территории Украины вегетационный период (период 
года, когда возможен рост и развитие растений) начинается в конце марта, а закан- 
чивается в конце октября, и наиболее присущими являются летние месяцы, то для 
анализа температуры целесообразно выбрать среднюю температуру в летний период 
(например, среднюю температуру воздуха в июле месяце) [8]. Для этого удобно исполь- 
зовать карту соответствующих изолиний температур, созданную с помощью средств 
МаршюЮ (линии, полилинии, полигоны). 

Получив карту изолиний средних температур в июле месяце для территории 
Украины, формируем карту лишь для анализируемых областей. Запрос языком МарВазс 
выглядит следующим образом: 
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Мар Егот 1етр_р1аруз, ВезВоипа_ НаОКг 

Зе Мар Гауег 1 Е4йаЫе Оп 

зеё тар гедга\ оп 

ЗеЕ Тагое! Оп 

зеё тар гедгауу о 

Зе Мар Огаег 2,1 

зе тар гедга\ оп 

ОБесё$ ПиегзесЕ шо Тагоеё Даа р127=р127, ЗОАТО 9=ЗОАТО_39, ОКГМате=ОКг Мате, 
КизМате=ВизМате, Еп®Мате=Епт?Мате, СешегОКгМате=СещегОкг Мате, 
СешегЗОАТО_О=СещегЗОАТО_Ч, Рори! 97=Рори197, ОтбапРори!97=ОгфапРору197, 
Кига[Рори!97=Вига!Рори! 97, Рору198=Рори198, ОтбапРори!198=ОгБапРорч198, 
Кига[Рори!98=Вига!Рори198, Тегпрегавхе=Тетрегайиге 


Показанная на рис. 3 цифровая карта отображает расположение изолиний 
средних температур воздуха в июле месяце для анализируемой территории Украины. 


— 
- 6 Чье 
№ та —— А 2 —% 
7 
\ . 
р. 


7 з\ ола 
> ? “ 
<. ы в. 


—19—— Изотермы средней температуры воздуха в июле месяце 


Рисунок 3 — Карта изотерм средней температуры воздуха в июле месяце 


Данные о типах лесов, температуре изолиний и другие, необходимые для анализа, 
параметры представляются в форме соответствующей таблицы. 


Формирование множества 
элементарных участков для анализа 


Для осуществления анализа климатических и физико-географических параметров 
для ареалов лесов на рассмотренной территории необходимо в первую очередь сфор- 
мировать элементарные участки, в пределах которых будет осуществляться такой ана- 
лиз. Для этого необходимо иметь множество элементарных объектов цифровой карты 
анализированной территории. За основу деления на такие объекты выбираем элементар- 
ные участки размером [ х [км (например, 5 х 5 км), ограниченные границами региона. 
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Пусть у, —-п-й элементарный участок (элементарный объект цифровой карты). 
Тогда множество всех элементарных объектов можно записать следующим образом: 


= ЬМ, (1) 


где №М- общее количество элементарных участков. 

Карта элементарных участков также может быть представлена в виде таблицы. 
Каждому элементарному участку соответствует строка цифровой карты с необходимыми 
данными и параметрами для расчетов. 

Реализовать процедуру деления карты региона можно, используя средства 
системы МарЮ по такой схеме: 

1) прибавить к открытой карте регионов объект «сетка» с соответствующими 
параметрами (размер клетки 5х5 км) и наложить ее поверх анализируемых 
областей; 

2) изменить свойства сетки и установить возможность редактирования; 

3) выделить сетку и прибавить объект, который изменяется, к выборке: 

Зе{ плар гедга\\ оЁР 

Зе Мар Гауег 1 ЕацаЫе Оп 
Зе пар гедга\и оп 

Зе! Тагое! Оп; 

4) карту региона переместить поверх сетки, выделить области для деления и 

удалить внешнюю часть, то есть часть сетки, которая очутилась за пределами региона: 
Зе{ плар гедга\и о 
Зе Мар Огаег 2,1 
Зе плар гедга\и оп 
ОБес$ Пиегзес{ шо Тагое аа Пезсирноп=Оезсгриоп, 
Со! Мате=Со! Мате, Во\и_ Мате=Ко\и Мате. 

Таким образом, получаем карту областей, которая разделена на элементарные 
участки. Такую карту можно использовать в качестве сетки для деления потенциальных 
лесов, а также для последующего анализа климатических и физико-географических 
параметров. 

Рассмотрим объекты, которые образуют цифровую карту лесов, и сформируем 


множество всех объектов леса. Пусть у —[-й объект цифровой карты лесов. Тогда 
совокупность таких объектов можно записать как множество: 


9, 1=ЬЕ)}, ©) 


где Г — общее количество объектов цифровой карты потенциальных лесов. 
Следовательно, для пространственной оценки расположения лесов рассматри- 
ваем два множества: 


М р ПУ — совокупность объектов, которые соответствуют лесу, из цифровой 
карты потенциальных лесов; 
|, ‚Ъп=ЬМ } — совокупность элементарных участков размера /х/ км. 


Теперь для описанных выше множеств применим операцию пересечения с расщепле- 
нием элементов. Результатом является множество 


о, 1 11; п Гм} о, 1 рп =ЬМ}, (3) 


(0 


п 


элемент у, определяет участок, который соответствует объекту вида потенциального 


леса в каждом элементарном участке. 


«Штучний 1нтелект» 32009 143 


Пасичнык О.Б. 


ЗП 


Таким образом, получаем карту потенциальных лесов, каждый объект которой 
разделен частью элементарного участка, в который он попал. Совокупность объек- 
тов на карте обозначим М, . 


Алгоритм анализа показателей температуры 
и абсолютной высоты 


Показатели температуры воздуха в июле месяце для анализируемого региона, а 
также абсолютной высоты, поданы в отдельных таблицах для каждого объекта в 
частности (изолинии для температуры и высоты). Но для проведения последующего 
анализа необходимо найти климатические данные для каждого элементарного участка. 

Для этого необходимо выполнить следующее: 

1) на открытые карты изолиний температур в летний период и абсолютных 
высот наложить сформированную карту элементарных участков; 

2) добавить к таблице элементарных участков две новые колонки, в которых 
дальше будут записаны соответствующие значения высоты и температуры; 

3) в таблице элементарных участков анализируется каждый из них, и если 
через него проходит изолиния высоты или температуры, то записываем фиксированное 
для элементарного участка значение параметра в соответствующую строку таблицы; 

4) в клетки, через которые не проходят изолинии, записываем значение 0; 

5) определяем наименьшее значение координаты У для карты элементарных 
участков; 

6) анализируем постепенно все элементарные участки, координаты У которых 
одинаковые, то есть рассматриваем каждую строку карты, двигаясь слева направо и 
сверху вниз; 

7) находим два последовательных значения одного климатического параметра, 
между которыми есть несколько незаполненных нулевых значений; 

8) равномерным делением заполняем все промежуточные нулевые клетки, 
причем если зафиксированные значения последовательных параметров растут, то и 
промежуточные значения будут постепенно расти, и наоборот; 

9) для элементарных участков, которые остались незаполненными, прорабаты- 
ваем такую же процедуру, но двигаемся сверху вниз, выбирая столбцы с одинаковы- 
ми координатами Х. 

Выполняя такую процедуру для каждого климатического параметра, получим 
карту элементарных участков с соответствующими значениями средней температуры и 
абсолютной высоты в каждом элементарном участке. Но для последующего 
исследования нужны климатические параметры не только для элементарного участка, но 
и для каждого вида леса. 

Для этого необходимо выполнить такие шаги: 

1) на открытую карту потенциальных лесов наложить карту элементарных участ- 
ков (с известными значениями параметров в каждой из них), в результате чего полу- 
чим карту потенциальных лесов, разделенную на элементарные участки; при этом 


получаем М№ элементарных объектов, каждый из которых помечаем у,, где 


п=ЬМ№М,1=1Г, ; в каждом элементарном участке имеем Г, объектов леса; 


2) в таблице потенциальных лесов, разделенных на элементарные участки, по 
очереди рассматриваем каждый элементарный участок и все объекты леса, которые 
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к ней принадлежат; сначала рассматриваем каждый объект леса по очереди в первом 
(п =1) элементарном участке и т.д.: у, п=1,1=1,[, ; 


3) записываем среднюю температуру и абсолютную высоту элементарного участка 


для каждого элементарного объекта леса 1=1,[, в рассмотренном элементарном 
участке у, п=1: 


АЗ СоГлллп Еогезё Огуае (Н1еЪ) Егот КИта Ра Зе! То Н12ВАБ$ 
МПеге пиегзесв$ Оупапис 

АЗ Со[лтлп Рогезё Пгуе (Тетрегаехге) Егот КИта Да Зе! То 
ЗиттегТетрегаиге У/Веге пиегзес{5 упатс 

Получаем карту потенциальных лесов, разбитую на элементарные участки, и 
для каждого объекта леса в каждом элементарном участке имеем значение температуры 
в летний период и абсолютной высоты. 

В результате выполнения таких действий получаем множество объектов цифровой 
карты, каждый из которых соответствует элементарному участку леса определенного 
типа, отнесенному к соответствующей клетке сетки размером 5х5 км. Общее число 
объектов цифровой карты составляет 


Ме У В (4) 


[] Уедуре Уедзреез Оезспр@ол |Агеа Тетрегаиге | Нов 4$ 

[9] 1 300 | Карпатские горные и предгорные леса | Ель, пихта, бук, дуб 8-38 25,094410 18,00 400,00 
"ЕЯ 1 300 | Карпатские горные и предгорные леса | Ель, пихта, бук, дуб В0103 25,203893 18,00 1 000,00 
о 
О 
О 
О Ель, пихта, бук, дуб в6105 23.549521 18.00 1 000,00 
| Ель, пихта, бук, дуб 25,226415 18,00 
| 


1 200 | Широколистные леса Бук, дуб, раб, клен, липа |В648 23,972476 1814 200,00 
1 200 | Широколистные леса Бук, дуб, раб, клен, липа |ВС47 23,949666 1816 200,00 


м 

ОГ зао [весоменьелья [Ба пабиленлиь 808 — | 8] 97| 00 
ОГ 92 [шяокомтные ля [БА пябимилиь [608 | м] #5 м 
ОГ 92 [швоомоные ля — [БА пабиинлиь 604 | 2298] _ #5 20 


О 
41 


1 200 | Широколистные леса 


Бук, дуб, раб, клен, липа |ВР52 18,755335 


1118336 18,08 


225.00 


Бук, дуб, раб, клен, липа |В653 


Рисунок 4 — Фрагмент таблицы с информацией о климатических и физико- 
географических параметрах для каждого объекта леса в элементарном участке 


Для анализируемых областей Украины количество объектов цифровой карты 
составляет 


№ 
№ =» 2; =1101. (5) 

п=1 
Фрагмент таблицы, с информацией о климатических и физико-географических 
параметрах для каждого объекта леса в элементарном участке, а также о других 

необходимых данных, представлен на рис. 4. 

Для анализа влияния климатических и физико-географических параметров на 
ареалы лесов разработана геоинформационная технология, которая базируется на создан- 
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ных программных модулях с использованием языка СЯ и системе геораспределенной 
базы данных. Результатом работы созданного программного обеспечения является 
обобщенная информация о наиболее благоприятных условиях прорастания основных 
видов лесов. 


Выводы 


Разработанная математическая модель, отображающая влияние физико-географиче- 
ских и климатических параметров на ареалы лесных пород, дает возможность получать в 
обобщенной форме информацию о наиболее благоприятных условиях прорастания ос- 
новных видов лесов. Созданный инструментарий демонстрирует возможность реализа- 
ции геоинформационной технологии, которая в сочетании со средствами моделирования 
сценариев изменения климата даст возможность анализировать и прогнозировать влияние 
изменения климата на ареалы лесов, давать рекомендации по ведению лесохозяйственной 
деятельности. 
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О.Б. Паачник 

ГеонформашИна система аналзу клматичних 1 ф1зико-географ1чних параметра для досл\дження 
ареалв лених порд 

У статти розглядаеться гео1нформащйна система аналзу кл1матичних (середня температура) та флзико- 
географучних (абсолютна висота) параметрав для дослдження ареалв л1сних порд. Запропонований шдхд 
до анал!зу параметрав, який дозволяе моделювати сценарий зм1ни кл!мату та прогнозувати вплив змни 
кл!мату на ареали лев. 


О.В. Раясйпук 
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